Edge states and enhanced spin-orbit interaction at graphene/graphane interfaces 
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We study interfaces between graphene and graphane. If the interface is oriented along a zigzag 
direction, edge states are found which exhibit a strong amplification of effects related to the spin- 
orbit interaction. The enhanced spin splitting of the edge states allows a conversion between valley 
polarization and spin polarization at temperatures near one Kelvin. We show that these edge states 
give rise to quantum spin and/or valley Hall effects. 
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Two-dimensional electronic systems in which the im- 
portant physics is taking place solely at the edges have 
received considerable interest during the last decades, 
starting with the discovery of the quantum Hall effect 
(QHE)R In a strong magnetic field, gapless transport 
channels exist only at the edges of the system and these 
give rise to a remarkably well quantized conductance. For 
a long time, the QHE was the only known effect of this 
type. Recently, however gap less modes have also been 
predicted^ and observecP^ on the surface of topologi- 
cal insulators without the need of magnetic fields. The 
first syste m in w hich these states have been investigated 
is graphen d 2 * 6 * 7 [ Here, spin-orbit interaction (SOI) plays 
the role of the magnetic field in that it opens a bulk 
energy gap. Spin-polarized transport is confined to the 
graphene edges which are oriented along the zigzag direc- 
tion. It turned out, however, that the SOI in graphene is 
rather wealP. Furthermore, it is difficult to actually con- 
struct clean structural edges of graphene, although some 
experimental progress has been reported^. 

Recently, a new possibility for the construction of 
very clean graphene terminations has been proposed^ 
rather than structurally cutting graphene in order to 
create nanostructures, it has been proposed to hydro- 
genate it locally. Thereby, the 7r-band is re moved wher- 
ever graphene is transformed into graphan d 11 * 12 * and an 
effective edge is created for the 7r-electrons. By this tech- 
nique of local hydrogenation, not only ribbons with high 
quality edges of the usual zigzag type (we call them ex- 
edges, henceforth) are possible; also bearded edges (/3) 
are within reach (see Fig. [TJ. 

We study such graphene/graphane (GG) interfaces of 
a- and /3-type and find edge states which are exponen- 
tially localized at the interface. However, instead of be- 
ing energetically nearly flat, like in the case of structural 
graphene edges, the edge states at GG interfaces have a 
considerable dispersion, the amplitude of which is largely 
determined by the energy of the hydrogen ls-orbital. 
These edge states enable a quantum spin-valley Hall ef- 
fect (QSVHE). Moreover, the SOI of these edge states is 
strongly enhanced due to the proximity of graphane. As 
an application we propose a device which is capable of 
converting valley polarizations into spin polarizations. 

Edge states. Henceforth we assume that all edges and 
interfaces are oriented along a zigzag direction. Due to 



the two-atomic unit cell of a hexagonal lattice, there are 
in principle two different types of boundary conditions 
- cv-type and /3-type (see Fig. [I]). The reason why a 
/3-edge is seldom considered is that it is highly unsta- 
ble against structural recombinations^. Indeed, we are 
not aware of any experimental observation of a stable fl- 
edge in graphene^. However, if we do not require the 
graphene flake to be cut in order to generate a /3-edge, 
but only require the hydrogen deposition on graphene to 
start on a certain sublattice while the hexagonal struc- 
ture remains intact, it makes not much of a difference 
whether a- or /3-interfaces are to be created. 
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FIG. 1: (color-online) a- and /3-interfaces in a 
graphene/graphane heterostructure. The large (red) 
spheres represent the carbon atoms and the small (blue) 
spheres the hydrogen atoms. The brighter spheres correspond 
to the graphene region while the darker spheres correspond 
to graphane. 

Before we turn to a more elaborate modeling of GG 
interfaces and the related edge states which respects the 
7r-band and the a-band, we would like to start with a 
qualitative discussion on the basis of a very simplified 
model which, as it turns out, captures most of the im- 
portant concepts. We begin with a nearest neighbor 
tight-binding model for the 7r-electrons in graphene Hq = 
*E(rr') c r c r' where t ~ — 3eV, (r,r') runs over nearest 
neighbors of a hexagonal lattice r = niai + n 2 a 2 + sRi 
(ni,»i2 € Z, s — 0, 1) with ai,a2 the Bravais lattice vec- 
tors and Ri the vector connecting the A and B sublat- 
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tices. c r is a 7r-electron annihilation operator. We want 
to describe geometries which are lattice-translationally 
invariant along a zigzag direction, which we choose to 
be parallel to ai (henceforth we call this direction the 
x-direction) . Thus, it is convenient to transform to elec- 



tron operators d n 



k,s 



where 



N x is the number of unit cells in the ribbon along the 
x-direction, so that the ID Brillouin zone is [0,27r]. In 
terms of the d-operators, the Hamiltonian reads 



projected Hamiltonian reads 



w,° 



o 

t' 



(4) 



In leading order perturbation theory in and ^, this 
projection yields a low-energy edge state \ip a (k)) ~ 
H \^{k)) - iA/£ \H) + 0(e H ) \C) with energy 



H G = t ^2 d n,k,A d n,k,B 
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where n labels the position along the a 2 direction and wnere wc h ave dropped terms of order fft and Note 



Uk = 1 ■ 



A structural a-edge is created by removing 



all terms with operators corresponding to positions n < 
and the terms with operators corresponding to the A 
sublattice site of n = 0. For a (3-edge only the n < 
terms need to be removed. A Hamiltonian truncated in 
such a way gives rise to exponentially localized, zero- 
energy edge states^ 



^ /P {k))=Nl 



with the normalization Af^ 



" J n,k,BIA 



(1 - K| ±2 )5, the local- 
ization length = =p [In 1 and some unimportant 
phase 4> = ar g( — u fc 1 )- Obviously, the a-edge state only 
exists for fee [^f, ^f] while the remaining k-space sup- 
ports the /3-edge state. Exactly at the projection of the 
valleys K and K' to the one-dimensional Brillouin zone, 
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i.e. ^ and K' 



4^ , the localization lengths of both 
2 iJfc-K'l- 1 , re- 
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edge states diverge as — K| 1 or 
spectively. For a ribbon with a- and /3-edges on opposite 
sides, a zero energy mode is present throughout the whole 
Brillouin zone. The corresponding wave function, how- 
ever, 'jumps' from one edge to the other as K or K' is 
crossed. 

Now, instead of a structural graphene edge we want to 
describe a GG interface (we discuss only the a-interface, 
the /3-interface being analogous, see Appendix A). The 
essential features are well described within a simplified 
model which includes only one -C-H group at each 'tooth' 
of the zigzag edge, instead of the whole graphane lattice. 
We further neglect all er-orbitals of the carbon atoms. 
The Hamiltonian of this -C-H group reads 



Hr = td 



'k,c d 0,k,B + t'dl H dk,c + h.c. + £nd\ H d k ^ H (3) 

and must be added to the a-truncated Hq, where d k ,c 
and d k ji are the annihilation operators of the carbon 
7r-orbital and the hydrogen s-orbital, respectively. Typ- 
ically, €h — —OAeV is much smaller than the other two 
energy scales t ~ —3eV and t' ~ — b&eV (see Appendix 
B). Since we are only interested in how the -C-H group 
affects the edge state, we project Hq + Hi onto the sub- 
space {m(k)),\C) = 4 >c |0),|#) - d\ H |0)}. The 



that the amplitude of the state \ip a (k)} is of order en at 
the C atom of the -C-H group. As a result, the inclusion 
of more and more graphane rows, i.e. -C-H groups, only 
leads to terms in the energy which are higher order in 
£h- Thus, this simplified model is expected to describe 
the energy dispersion correctly as long as no other states 
which we have neglected here come too close in energy 
Eq. ([5]) roughly describes a parabola around k = n which 
crosses zero at K and K'. Therefore, the edge state at 
a GG interface has - in contrast to the edge state at a 
structural zigzag edge - a finite velocity, the magnitude of 
which is largely determined by the energy of the hydrogen 
ls-orbital e# relative to the 7r-orbital energy in graphene 
(set to zero here). Note that also other effects which are 
not included in our simplified model, e.g. next-nearest 
neighbor hoppings^ or local electrostatic gates at the 
edges^, can increase the bandwidth of the edge state. 
For typical tight-binding parameters we find edge state 
velocities near K or K' of 10 4 — 10 5 ™ (see Appendix A) . 

The bandwidth of the edge state is crucial for the ques- 
tion of magnetic ordering at the edges. While conven- 
tional graphene ed ges ar e believed to be spin-polarized 
in their ground stat e 1 18 1 19 - 1 , the edge states considered here 
have a larger bandwidth which helps to suppress the mag- 
netic phase transition. The bandwidth can be tuned fur- 
ther by gatesl 121221 . 

Fig. [2] compares the energy dispersion of the edge 
states calculated from the extended tight-binding model, 
as introduced below, with Eq. d5| and its analog for 
the /3-edge state (see Appendix A). While the simpli- 
fied model for the a-edge state agrees well with the ex- 
tended model, the simplified /3-edge state dispersion pre- 
dicts a too small bandwidth. However, the main features, 
namely the different directions of motion of electrons in 
the valley K (K') at different interfaces, are described 
properly by the simplified model. 

Edge state transport. A ribbon with combined a- 
and /3-interfaces (a/3-ribbon) with finite width W of the 
graphene region gives rise to a quantum valley Hall ef- 
fect (QVHE}2I1 in which K (K') valley electrons move 
left (right) at the a-interface and right (left) at the (3- 
interface, if the Fermi energy cf is tuned as indicated in 
Fig. [2] For wide enough ribbons (£fc F <C W/2), the over- 
lap between edge states moving in different directions is 
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FIG. 2: (color-online) The band structure of a graphane- 
terminated a/3-ribbon with W = lOnm. The dashed hori- 
zontal line is a typical Fermi energy €f for the QVHE regime 
chosen such that £fc F <C W/2. The gray bands correspond 
to exponentially localized states at the outer graphane edge 
which are not important here. The full curves are results 
of numerical calculations within the extended tight-binding 
model. The dashed (blue) curves show the effective model 
dispersion (|5|. The inset shows the conductance G of the 
graphene region as a function of ep. 

small so that backscattering is forbidden if valley scatter- 
ing can be neglected. Since valley scattering requires the 
transfer of a large crystal momentum, only atomic scale 
disorder would destroy the QVHE. 

The detection of this QVHE by a transport measure- 
ment could be a first experimental step to reveal the GG 
interface quality. In a clean a/3-ribbon the conductance 
is =#-n, where n is the number of transport modes inter- 
sected by e_p (see inset of Fig. [2]). 

For wide (W ^> 50nm; see Appendix A) a/3-ribbons at 
temperatures below the spin-orbit gap in bulk graphene 
(~ 10mK 8 ) a new quantized Hall effect arises: the 
QSVHE. Here, not only the spin or the valley determines 
the direction of motion at the edges, like in the QSHEtPor 
in the QVHE, respectively, but both, spin and valley are 
responsible for the selection of the edge and the direction 
of motion, if cf is tuned into the bulk SOI- induced gap in 
graphene. For instance, at an a-interface only electrons 
in valley K and with spin pointing into a certain direc- 
tion (we define this direction as f) are allowed to move 
left; the electron in valley K' and with spin pointing in 
the opposite direction (J,) moves right. The other two 
combinations of valley and spin are not allowed. This 
increases the stability of the QSVHE compared to the 
pure QSHE or QVHE because backscattering requires si- 
multaneous valley and spin scattering. In principle the 
QSVHE already allows spin-valley conversion, but only 
at very low temperatures (~10mK). 

Enhanced spin-orbit interaction. It is well known that 
the smallness of the SOI in graphene is rooted in the 
lightness of carbon on one hand and in the symmetry- 
induced decoupling of the ir- and er-bands, on the other. 
The latter issue can be overcome, however, by twisting 
the lattice. While the carbon atoms of the A and B sub- 



lattices are coplanar in graphene, they are 'pushed' out of 
the plane in graphane due to the rehybridization sp2 — > 
sp 3 which is caused by the presence of the additional 
hydrogen atoms^. Thus, the 7r-band is locally coupled 
to the (T-band at the interface. We therefore expect an 
enhanced SOI near GG interfaces. 

In order to substantiate this expectation, we use an 
extended tight-binding model of GG heterostructures 
which respects the 7r-band as well as the cr-band of 
the carbon lattice (see Appendices). It is based on an 
environment-dependent tight-binding model for general 
hydrocarbons^! and takes the atomic carbon SOI into 
accounlP. We find that the GG edge states exhibit a 
much larger spin splitting than conventional edge states 
(see Fig. [3|. This can be understood for the a-edge state 
(for the /3-edge state, see Appendix F) by projecting the 
on-site SOI Hamiltonian 

HsO=^AY / E ^ P < P ^T<r>Cr, Pp ,r', (6) 

where c r , p T annihilates an electron at site r in orbital 
p^ with spin r, are the Pauli matrices for the electron 
spin, fi, v, p = x,y,z, and e IJ,vp is the Levi-Civita ten- 
sor, onto the two-dimensional spin-dcgcncratc subspacc 
of the a-edge state \ip a (k); r) obtained by diagonalizing 
the extended tight-binding Hamiltonian without SOI. In 
a conventional edge state \ipQ (k); t) only the 7r-orbitals 
are occupied so that (ipg (fc); T\Hso\4>]i{k)\ t') = and 
the SOI becomes a higher order effedP. At a GG inter- 
face, however, the edge state acquires contributions from 
the cr-orbitals because of the out-of-plane twisting of the 
carbon atoms which leads to a mixing of the it- and the 
CT-band. As a result, the SOI becomes a first order effect. 
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FIG. 3: (color-online) SOI induced spin splitting of the edge 
states. The solid (blue) curve shows the absolute value of 
the spin splitting Aeso (k) of the edge states calculated from 
the extended tight-binding model which includes the bulk 
graphene SOI as well as the interface SOI. The dashed (red) 
line shows the spin splitting calculated from the approximate 
Hamiltonian j9j|. The solid curve has been calculated for a 
lOnm wide a/3-ribbon. The spin splitting at K,K' which is 
larger than the spin-orbit splitting in bulk graphene is due to 
finite-size effects. 
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The projected Hamiltonian reads 

H fo a = J2L ^ e L,r r rr'(fc)efc, a ,r', (7) 

TT' ~3~ 

with efc, a ,r the a-edge state annihilation operators and 

T TT ,(k) = {r{k)-r\H so \i> a {k)-T') (8) 
~ (fc-7r)(A^) 2 [A^ + A^] rT ,, (9) 

where A^,Ai are constants describing the Rashba and 
intrinsic parts of the effective SOL The Dresselhaus term 
(oc a x ) vanishes because of the mirror symmetry x — > —x 
of the interfaces. The additional factor (AT k ) 2 accounts 
for the fact that the effective SOI, generated by a GG 
interface, must be proportional to the amplitude of the 
state at the interface. The parameters A^ = — 0.16meV 
and Af = — 0.05meV are obtained from a fit to the 
numerical results (see Appendix F). Because of time- 
reversal invariance the spin splitting is exactly zero at 
k = 0,±7r. These are also the points where the spin 
direction of the energetically higher edge state changes 
abruptly. Note that for wide ribbons (W > lOOnm), the 
spin splitting at K,K' is essentially given by the bulk SOI 
because M k ^ vanishes for k -^K,K'. As a result, the SOI 
enhancement at the interfaces does not increase the crit- 
ical temperature at which the QSVHE can be observed. 



squeezing the wide edge states into a narrow ribbon, the 
spin-orbit splitting is enhanced at K and K'. 

As required, the energetically higher band at K has the 
opposite spin-direction than the corresponding band at 
K\ Thus, by tuning tf as indicated in Fig. [4] a spin- 
valley filter is realized: left- and right-movers exist for 
each spin and valley but the transmission of this struc- 
ture is only non-zero for (K,|) or for (K',|) and zero 
for the other two combinations of valley and spin. This 
means that each spin-up electron which initially consists 
of K and K' components (valley-unpolarized) will be in 
a pure valley state (K here) after it has passed the nar- 
row a/3-ribbon in x-direction, while a valley-unpolarized 
spin-down electron will be in a pure K' state after the 
passage. 

In conclusion, we have demonstrated that graphane 
can be utilized to create an effective edge for the 7r-band 
in graphene and that the SOI of the resulting edge states 
is strongly enhanced, compared to pure graphene. Fur- 
thermore, we have shown that, due to enhanced SOI at 
graphene/graphane interfaces, binary information which 
is encoded in the spin of an electron can be transferred 
to a valley-based encoding. 

We acknowledge useful discussions with B. Braunecker, 
I. Martin, B. Trauzettel, and A. Yacoby. This work 
has been supported by the Swiss NF and the NCCR 
Nanoscience Basel. 
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FIG. 4: Energy dispersion of the edge states of a 2nm wide 
a/3- ribbon near the K-valley. The two bands have opposite 
spin direction. The dashed line indicates a typical 6f- The 
inset shows a larger region of the band structure. Grey lines 
see caption of Fig. [2] 

However, there is another possibility to use the inter- 
facial enhancement of the SOI for high temperature spin- 
valley conversion. The spin splitting at K and K' decays 
like W~ l . Thus, one should use narrow ribbons in or- 
der to exploit the SOI enhancement at GG interfaces. In 
Fig. |4j the band structure of a narrow GG heterostruc- 
ture is shown. At the maxima of the edge mode the 
spin-orbit splitting is about 78 fieV ~ 0.9K. Note that 
although these states are derived from edge states, they 
do not actually appear like edge states in this case be- 
cause £fc 3> W (£fc diverges near K and K'). Thus, by 



APPENDIX A: SIMPLIFIED MODEL FOR THE 
EDGE STATE 

The a- and /3-truncated Hamiltonians for the struc- 
tural graphene zigzag edges read 

oo 



H G = [ d l.k,A d n,k,B + dl k A U k d n - ltktB \(Al) 

k n—1 
( oo 

H G = \ [ d l,k,A d n,k,B +dl !k>A U k d n - 



-l,k,B 



k \ n—1 

- d i,k,A d n,k,B^ 



(A2) 



with Uk = 1 + e ik and the simplified model Hamiltoni- 
ans for the a- and /3-interfaces between graphene and 
graphane read 

Hf = tdl c d 0ik ,B+t'dl H d k>c +h.c.+e H dl H d k:H (A3) 



#f = tu k d\ c d . k ,B + t'd\ H d ktC + h.c. + e H d\ H d k . H , 

(A4) 

where we have defined 

d k , H - N-?Y, e ~ iknic nuH (A5) 
til 

d k ,c = N-iJ^e'^CnuC- (A6) 
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The operators c nit c annihilate an electron in the tt- 
orbital of the carbon atom next to the edge, as shown 
in Fig. [5] and the operators c nit n annihilate an electron 
in the ls-orbital of the corresponding hydrogen atom. 




a-interface 



.-interface 



FIG. 5: Definitions in the simplified interface model. 

The green dots represent the carbon atoms of the -C-H groups 
and the blue dots represent the hydrogen atoms. The black 
numbers indicate the n\ coordinate of the carbon atoms in 
the graphene region. 



The projection of H^J^ 

|i) = 4, c l | 



Hf/P onto the subspaces 



{\rj"{k) 

matrix form 



. |2> = 4,jtI°>} 



reads 



m+Hf. 



proj. 



proj. 





t 

tu k M^ 



f 



t' 



(A7) 



(A8) 



Because |e#| \t\, \t'\, we are only interested in results 
to order eh . We hnd a low energy edge state for both a- 
and /3-interfaces 



t 2 {M^f 



t 2 (^k) 2 



t' 2 



e H + 0{e 2 H ), ke 



2n 4tt 



J* - 



* 2 KI 2 (0 2 



t 2 



Uk\ 



2 (A/-f) 2 



t' 2 



e H + 0(e 2 H ), k€ 



(A9) 
2?r 2tt~ 

'T' T 

(A10) 



Near the boundary of the k-space domains of the edge 
states, K= ^ and K'= 4^, the dispersion is approxi- 
mately linear and we can write for both states [a/ (3) in 
both valleys (K,K ! ) 

e a ^ = V3^\Sk\ + (D(Sk 2 ), (All) 



where 5k — k—K or 8k = k— K'. From this we can esti- 
mate the typical velocity of the edge states near K and 
K' to 



\v a/ P\ ~ ^3 



t 2 \e H \ a V3 
t 12 h 



6.8 • 10 4 



(A12) 



which is about one order of magnitude smaller than the 
Fermi velocity in graphene. 

Next, we would like to estimate the typical localization 
length of the edge states in the QSVHE regime. If we take 
the SOI into account, a spin gap A ~ 2/ieV opens up at K 
and K' for wide ribbons (W > lOOnm). This spin gap is 
of the order of the spin-orbit splitting in bulk graphene 
because the edge states become completely delocalized 
over the graphene region at K and K'. If the Fermi energy 
is of the order of the spin gap, the four Fermi momenta 
can be estimated by 



k x j 2 ~ K ± A 



K ± 3.4 • 10 7r. 



(A13) 



u3/ 4 ~ If' 



f'2 



K ± A 



V3t 2 e H 



K' ± 3.4 • 10" 3 tt. (A14) 



From this, the corresponding localization lengths for all 
interface states in the QSVHE regime are 



C P =2«o|lnK F ||- 1 ~22nm, (A15) 
with do — 1.4A the nearest neighbor C-C distance. 



APPENDIX B: TIGHT BINDING PARAMETERS 

We want to model the graphene/graphane heterostruc- 
tures by a nearest-neighbor tight-binding model which 
takes into account the 7r-band and the cr-band of the 
hexagonal carbon backbone and the ls-orbitals of the hy- 
drogen atoms attached to each C atom in the graphane 
region. We neglect the non-orthogonalities of orbitals 
on different sites. First of all, we need the bare tight- 
binding hopping integrals between the oriented carbon 
orbitals 2s, 2p x , 2pj,, 2p z , i.e. the matrix elements of the 
Hamiltonian H 



V ss = (2s;r |ff|2s;ri) 

V sp = (2p z ;r \H\2s;r 1 ) 

V; p = -(2 Pz ;r \H\2 Pz ;r 1 

V; p = (2 Px ;r \H\2 Px ;r 1 ) 



(Bl) 
(B2) 
(B3) 
(B4) 

and the hopping integrals between carbon orbitals and 
the hydrogen ls-orbital 

W ss = (ls;r |if|2 S ;r 2 ) (B5) 
W sp = (2p z ;r \H\ls;r 2 ) , (B6) 

where r = (0,0, 0) T , r x = a o (0,0,l) T and r 2 = 
&o(0,0,l) T with ao the nearest-neighbor C-C distance 
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and bo the C-H distance. The kets 1 2s; r) , \2pi; r) rep- 
resent the 2s- and the 2pi-orbitals (i = x, y, z) of the 
carbon atoms at r. The ket |ls;r) represents the hydro- 
gen ls-orbital. 

The H-H distance in the chair conformation of 
graphane is large enough to neglect the direct hopping 
between hydrogen orbitals. All nearest-neighbor hopping 
integrals between the relevant orbitals in graphane can 
be reduced to the parameters defined in Eqs. ( B1|B6 I, 
as we will show subsequently. 

In general, these bare hopping parameters are 
environment-dependeni)- 22 l 26 l However, there is a consid- 
erable variance of hopping parameters in the literature. 
For instance, the very often used parameters in Ref. HJJ 
deviate considerably from the ones extracted from Ref. 
|2"B1 (see Tab. [i]). In this case, the deviation is probably 
due to the different scopes of Refs. US] and [571 namely 
the crystal structure and the electronic properties, re- 
spectively. Therefore, we believe that the carbon orbital 
hopping parameters of Ref. [57] are more suitable to our 
needs. Furthermore, as Reich et al. argue^ for the 7r- 
band, for a quantitative description one would have to 
include up to 3rd neighbor hoppings and the correspond- 
ing non-orthogonalities. These parameters arc, however, 
to our best knowledge, not available for the cr-band in 
graphene. Since we only aim at a description on a quali- 
tative level, we refrain from more quantitative first prin- 
ciples calculation. We rather check that our results are 
robust with respect to variations in the parameters. 





Ref. [26] 


Ref. [27j (used here) 


Vss 


-4.6 eV 


-6.8 eV 


y PP 


-2.4 eV 


-3.0cV 


v pp 


-7.7eV 


-5.0eV 




-5.7eV 


-5.6eV 



TABLE I: Two sets of nearest-neighbor tight-binding param- 
eters from the literature. We finally use the second column of 
parameters in the calculation. 

Furthermore, we make the simplifying assumption that 
the bare hopping matrix elements arc the same in the 
graphene and in the graphane region. This does not 
mean, of course, that the final Hamiltonian in the carbon 
subspace is equal for graphene and graphane. The differ- 
ent relative alignment of the A- and the B-sublattice is 
taken into account later. 

The parameters for the C-H bond in the graphane re- 
gion (W ss and W sp ) are calculated along the lines of Ref. 
H3 with a C-H bond length of l.lA^. They are 



W s . 



-5.4eV, 



-5.8eV. 



(B7) 



Finally, the 2p-orbital energy of carbon defines the zero 
of the energy. The relative energy of the carbon 2s-orbital 
is e s c = —8.7 eWA The orbital energy of the hydrogen ls- 
orbital can be roughly estimated from Refs. 22 26. We 



assume that the orbital energies of the carbon orbitals 
is elevated due to the nearby hydrogen atom about the 
same amount as due to the other carbon atoms. In this 
work we use e s H = — 0.4eV for the hydrogen ls-orbital. 
The results and conclusions do not depend critically on 
this value. 



APPENDIX C: ORBITAL RIBBON 
HAMILTONIAN 

We first want to derive the Hamiltonian, describing 
bulk graphane and bulk graphene, in order to check the 
resulting band structure against ab initio band struc- 
tures in the literature^!. For doing this, we need to 
translate the bare hopping integrals ( B1|B7 1 to the hop- 
pings on the lattice under consideration. We define 

the honeycomb lattice vectors ai = clq (y3, 0, 0) and 



§, 0^ and the nearest neighbor vectors 



R 2 



ao(0, l,z) T 

V3 1 



«o 



2 ' 2' Z 



'3 1 



a V~2' 2 . 



(CI) 
(C2) 

(C3) 



where ao is the C-C distance of flat graphene and zclq 
is the separation of the A- and B-sublattice planes. We 
chose the parameter z = zq ~ 0.42 in this work. At 
each carbon site, we define the system of the 2s- and 2p- 
orbitals in the following way: the s-orbital is spherically 
symmetric, so that the atomic alignment is irrelevant, 
the p-orbitals can be characterized by a vector, pointing 
into the direction of the positive part of the orbital wave 
function. The wave functions of the second carbon shell 
are approximately 



fe(r) = / 2s (|r|) 
^2 Pi (r) = /2p(|r|)rj, i 



x,y,z. 



(C4) 
(C5) 



Thus, we choose the alignment of the three 2p-orbitals 
according to the same coordinate system which defines 
a 1; a 2 . The set of p-orbitals transforms as a vector. This 
is useful if we want to express a p-orbital pointing, say, 
into the (1,1,0) direction. The wave function of this or- 
bital can be written in terms of the basis functions ( C5 1 



^2 P , (i,i,o) (r) = [fp2p x ( r ) + V>2 Ps ,(r)] 



(C6) 



The direction of an orbital is conveniently expressed by 
a vector p. 

We need to deal with three different types of C-C hop- 
pings. There is the hopping from an s-orbital to another 
s-orbital (ss) , the hopping from an s-orbital to a p-orbital 
(sp) and the hopping between p-orbitals (pp). For the 
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s-orbital 



orbital 



s-orbital 




p-orbital 







p -orbital 











FIG. 6: Relative orientation of the carbon s- and p-orbitals. 
rt is the bond vector, connecting the positions of the two 
atoms. 



subsequent discussion we introduce the bond vector rj, 
which connects the two atoms participating in the bond 
under consideration. For an sp-bond points always 
from the atom, carrying the s-orbital, to the atom on 
which the p-orbital is located. For ss-bonds and pp-bonds 
the direction of r;, plays no role. 

The ss- hopping does not depend on the direction of 
but only on the bond length. The bond length, however, 
is assumed constant over the whole lattice. In a homo- 
geneous structure (only graphene or graphane) , this as- 
sumption is fulfilled perfectly. In a heterostructure like 



the one we aim to describe there is a minor difference in 
the bond lengths in the different regions. We neglect this 
difference in the bond length because we believe that the 
dominant difference between graphene and graphane is 
the different symmetry of the lattice. 

For sp-hopping we only know the hopping integral if 
the p-orbital is aligned along the bond vector. The hop- 
ping between an s-orbital and a p-orbital that is perpen- 
dicular to the bond vector is zero by symmetry. Thus, 
only the angle 9 between the p-orbital and the bond vec- 
tor is important (see Fig. [6]). The sp-hopping integral of 
such a bond is then 



t sp (6) — —V sp cos 9, cos 9 = 



r& p 

MIpI 



(C7) 



For pp-hopping, three angles are relevant: the two an- 
gles 9\,9 2 between the p-orbitals pi,P2 and the bond- 
vector Tb (see Fig. [6| and the angle p between the planes 
spanned by r^, pi and rj,, P2- The hopping integral of 
such a bond is 

tp p (9i,d2, p) = cos 9 1 cos 9 2 Vp p + sin 9\ sin 9 2 cos pVp p 

(C8) 

with 

cos9i — —- t — - — -, cos 9 2 — ] — n — r 



m\Pi\ 



\ r b\\P2\ 



and 



p = <[r b x p!,r b x p 2 ]. 



(CIO) 



With the above equations we can evaluate all hopping 
matrix elements we need to write down the Hamiltonian. 
We start with the Hamiltonian in the carbon subspace 
He- We denote the sublattice by s — A,B, the Bravais 
lattice vector R n = niai + ri2a2 and the orbitals by 
p = 0, 1,2,3 for the s-, p^-, Pj,-, p z -orbital, respectively. 
The real-space Hamiltonian reads 



He — ^ C li,fi,A [^,M',l c n:M',-B + */i,/i',2Cn-(0,l),ju',B + ^^':3 c n+(l,-l)^',s] + h.C. , 



(Gil) 



where c njAljS arc electron annihilation operators and j for j — 1, 2, 3 are the hopping integrals between orbitals p 
and p! calculated by Eqs. (C7),(C8l for a bond between a central A atom and its three neighboring B atoms. For a 



bulk system, this Hamiltonian can be transformed to k-space as usual 



Br 



EOw'W^'.fl + h - c - 



r 



(C12) 



where and 
/ W '(k) = Wu + W'.ae -4 * 3 + tw^-M (C13) d k ^ s 



\ " e -i(kin 1 +k 2 n 2 ) c 



(C14) 
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Note that ki , k^ are the components of the k-space co- 
ordinate w.r.t. the reciprocal vectors, which are non- 
orthogonal. In terms of k x and k y we have ki — \/ok x 
and &2 = \{V3k x + 3k y ). 

The graphcnc regions of the hcterostructures are fully 
described by He- In the graphane region, the hydrogen- 
related terms must be added to He- Because we only 
consider nearest-neighbor hopping and the chair confor- 
mation of graphane, direct inter-hydrogen hopping is not 
allowed. Since the hydrogen atoms are supposed to sit 
nicely on top (z > 0) of the B sites and below the A sites 



(z < 0), the only non-zero hopping integrals, allowed by 
symmetry, are the two between the hydrogen ls-orbital 
and the carbon 2s- and 2p 2 -orbitals, i.e. W ss and W sp , re- 
spectively. Because of this perfect alignment, no further 
transformation of the hopping integrals is needed. Only, 
because the 2p 2 -orbital points in positive z-direction at 
the A- as well as at the B-sites, we have to respect the dif- 
ferent sign of the sp-hopping between the carbon atoms 
on A-sites and their attached hydrogen atoms and those 
on B-sites. We find 



H-C — 



k,s 



+ (-!)' 



ik,3,s 



(C15) 



r 



The spinless bulk Hamiltonian H = He + Hh-c can 
be represented as a k-dependent (10 x 10)-dimensional 
matrix which is easily diagonalized. Fig. [7] shows the 
graphane band structure calculated from H. Compared 
to the first-principles band structure of Sofo et aZpSl, the 
band width of the lowest band is somewhat larger in 
our tight-binding model. This is typical for tight-binding 
models on honeycomb lattices with only nearest-neighbor 
hoppings and orthogonal orbital wave functions. The 
band gap is with E g ~ 4.9eV also larger than the one 
found in Ref. [551 However, Lebegue, et alW%, point out 
that the band gap is rather 5.4 eV. The band gap of 3.5 
eV, as found by Sofo, is, according to Lebegue et al, 



due to the use of the generalized gradient approximation 
which is known to predict erroneous energy gaps. 

For modeling graphane-terminated graphene nanorib- 
bons which are lattice-translationally invariant along 
the ai-direction (zigzag ribbons), we perform a partial 
Fourier transform of the carbon-electron operators, i.e. 



1 



\ ' 



-ikn . 



(C16) 



The ribbon Hamiltonian in the carbon subspace then 
reads 



He 



E 

k 1 n,fj,,fj, / 



fi.fi' ,l^n.k,fi' ,B 



e lk t„ 



-l,k,n',B\ 



(C17) 



It is important to note that the hopping integrals t^^tj 
are different in the graphane and graphene regions. At 
the interfaces we assume that the carbon atoms in the last 
graphane row have the sp 3 -like orbital configuration of 
graphane rather than the sp 2 configuration of graphane. 
There is some ambiguity in this choice. Different choices 
of the interface properties do not lead to qualitatively 
different results. Only the magnitude of the interface- 
induced spin-orbit splitting may be renormalized by fac- 
tors of order one. 



The hydrogen part of the Hamiltonian is easily added 
to the graphane regions, as explained above. 



APPENDIX D: SPIN-ORBIT HAMILTONIAN 

We restrict ourselves to the on-site spin-orbit interac- 
tion generated by the Hamiltonian 



Hso = 



Am 2 c 



=A 



r(W(r) xp)-<x, 



(Dl) 



where ^(r) is the (a priori unknown) potential of the car- 
bon ions, p = —ih(d x , d y , d z ) is the momentum operator 
of an electron, and cr — (a x : a v , <r z ) are the Pauli matri- 
ces for the electron spin. We are interested in the matrix 
elements of Hso m the subspace spanned by the car- 
bon 2s- and 2p-orbitals, i.e. the states |Vv.r) = Im> t ) 
it (g) |r) where [i = 0,1,2,3 labels the s-, p^-, p y -, p z - 
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FIG. 7: Bulk band structure of graphane calculated from our 
tight binding model (Eqs. \C12\ and (|C15||). 



orbital, respectively, and r =t, J, labels the spin. 
( M) r| J ff so | M / ! r')=A (HTOxpl.lf*')^- 

(D2) 

Now we perform an explicit calculation of the matrix 
elements and relate all non-zero integrals to each other 
by symmetry considerations. In real space, we can write 
(i,j,k = x,y,z) 

( M | [W(r) x p] 4 = -ie^J d 3 rV^(r)(9^)9 fc ^(r), 

(D3) 



where we used the wave functions defined in Eq. ( |C5| . 
From the symmetry of the orbital wave functions and the 
potential V(r) it is easy to see that 



(s\H so \Pi) = { Pl \H so \s) = 



A(p, 



[Wxpl 



Pk 



= ie ijk A, 



where e 1 ^ is the Levi-Civita tensor and 

' drr\f p (r)fV'(r). 



Am 2 c 



In second quantization the Hamiltonian reads 



(D4) 
(D5) 



(D6) 



(D7) 



TT 



where jj,,v,rj = x, y, z. The coupling constant A ~ 3meV 
is determined from the atomic spin-orbit interaction (see, 
e.g., Ref. ®). 



APPENDIX E: GRAPH ANE-TERMINATED 
a/3- RIBBONS 

As mentioned above, there is some ambiguity in the 
alignment of the interface bonds. There are two extreme 
cases. The first is shown in Fig. [TJ there, the bonds be- 
tween the graphene and the graphane region are aligned 



as if the first graphene atom was sp 3 -hybridized, i.e. they 
are tilted out of the xy-plane just as all other bonds in 
the graphane region. The other extreme case would be 
that the interface bonds are completely in-plane, as if 
they would correspond to the graphene region. It is not 
known to which extent the bonds are tilted out of the 
plane. We believe that the situation is more like shown 
in Fig. [T] because there should be a kind of a repulsive 
force between the hydrogen on the first graphane atom 
and the first graphene atom. Ab initio or structure cal- 
culations based on environment-dependent tight-binding 
models could yield a better estimate of the interface de- 
tails. 

In the following we quantify the tilting at the interface 
by a parameter zj. The fully tilted situation (as shown 
in Fig. [T]) is described by setting zj — z . The flat 
interface, on the other hand, is described by zj = 0. 




FIG. 8: Different band structures for zj = (red), Zj = zq/2 
(green) and zj = zo (blue). The graphene part of the a/3- 
ribbon is 21nm wide in this calculation and the graphane 
terminations are 2.1nm wide. The bulk states and the states 
at the outer graphane edges are not significantly affected by 
zj. The corresponding energy bands lie on top of each other 
for the three cases. 

The impact of the interface details on the dispersion 
of the GG-interface state is shown in Fig. [8] Obviously, 
primarily the edge states are affected by the interface de- 
tails while the bulk dispersion is largely invariant. Also 
it is observed that, while the a-edge state's dispersion al- 
ways resembles a parabola very closely, the /3-edge state 
changes the shape of its dispersion more heavily. How- 
ever, for this work only the dispersion near K and K' is 
important. There, changing zj mainly renormalizes the 
interface state velocity by a factor of order unity. 

In the numerical diagonalization of the a/3-ribbon 
Hamiltonian the dispersion of the /3-edge states are some- 
times (depending on the interface details) crossed by 
bands which are exponentially localized at the outer 
graphane edge (see Fig. |8|. Usually, this crossing hap- 
pens at |fc| < 0.27T. These states are spatially separated 
from the edge states at the interfaces and are thus not 
important here. If the graphane regions were infinitely 
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wide, these bands would not affect the /3-edge state at all. 
However, since the numerical calculations use a finite size 
ribbon, these graphane edge states have an exponentially 
small overlap with the /3-edge state wave function. Since 
we are generally interested in the band structure near K 
and K', this crossing is not important for our reasoning. 
Thus, in the plots of the numerical spin splitting and the 
spin-orbit strength (see below), we leave out the part of 
the Brillouin zone which is beyond the crossing. This is 
why the plots of some numerical results are restricted to 
the /c-range [0.27T, 1.87r]. 



APPENDIX F: EFFECTIVE SPIN-ORBIT 
HAMILTONIAN 

As explained above, the effective spin-orbit Hamilto- 
nian of the edge states at the GG interface is obtained by 
projecting Hso (see Eq. (D7l) to the restricted Hubert 
space of the edge states \ip a {k); r). This leads to an ef- 
fective SOI described by 





Zj =0 


Zj = 20 




-0.048 meV 


-0.16 meV 


A? 


0.026 mcV 


-0.050 mcV 


K 


0.44 meV 


0.13 meV 


Af 


0.057 meV 


-0.045 meV 



TABLE II: Single parameter fit of the spin-orbit interaction 
for the edge states at a- and /3-interfaces. 





zj=0 


Zj = 2 


a£ 


0.56 meV 


0.056 meV 


< 


-0.09 


0.99 




5- 10" 5 


-0.18 


Af 


0.18 meV 


0.09 meV 


4 


-0.40 


-0.89 




0.05 


0.14 



TABLE III: Three parameter fit of the spin-orbit interaction 
for the edge state at a /3-fnterface. 



T T r'(k) = ^ a {k)-T\H so \^ a {k); 



(Fl) 



which can be decomposed by means of spin Pauli ma- 
trices. It turns out that the function, multiplying <j x is 
zero while the coefficient functions of a v ^ z are finite. This 
makes sense because in addition to the intrinsic SOI of 
pure graphene, there is a z — > —z symmetry breaking 
at the interface which leads to a Rashba-like term k x a y . 
The x — > — x symmetry, however, is not broken by the 
interface so that the Dresselhaus-like term k x a x vanishes. 

The function T TT i(k) can be calculated numerically 
from the wave functions of the GG edge states. How- 
ever, it can well be approximated by an analytic form: 
the spin-orbit interaction at the GG interface is strong 
compared to the spin-orbit interaction in bulk graphene. 
Therefore, an effective model must respect the renormal- 
ization of the SOI by the amplitude of the interface state 
wave function directly at the interface. This is given by 
Afu (A/jf ) for a- (/?-) edge states in infinitely wide rib- 
bons. 

Moreover, because of time-reversal invariance, the fc- 
dependent coupling strength must be odd at k — 0, tt. 
Thus, we only take into account odd polynomials around 
these two points. As it turns out, the linear terms are 
sufficient to fit the numerics of the a-edge state. For the 
/3-edge state more terms are needed for a quantitative fit. 

The guessed form of the effective spin-orbit interaction 
of the a- and /3-edge states with only linear terms reads 



-pcff 

1 so 



a (k) = (k~n)(^) 2 [A a R ay + Afa z ] (F2) 



«//3 



and A"^, we fit 



In order to find the parameters A R 
these functional forms to the numerical results. Fig. [9 
shows various fits. In Table |H| the effective spin-orbit pa- 
rameters for the flat (zj = 0) and the tilted (zj — zq) 



interfaces are given. The a-edge state spin-orbit interac- 
tion is well reproduced by the effective form given in Eq. 
(F2| while the /3-edge state spin-orbit interaction is not. 
Therefore, we include higher powers in Eq. (F3 1 



-ieff 

so, p 



(k) 



M<) 2 



flC 1 t "fl 
+ Af(l + affc 2 + 6ffcV 



(F4) 



The parameters fit to the numerical results are given in 
Tab. [II] for the one parameter forms (|F2| and (LF3| and 
for the three parameter form (|F4| in Tab. Ill In order to 



compare the effective models with the numerical results, 
we introduce the Rashba (intrinsic) spin-orbit strength 
Tft(k) (Ti(k)) as the coefficient of the Pauli matrix a y 
(cr z ) in the expressions ((fTJ, ( |F2| , ( |F3| , and flF4| . Fig. 
[9] compares the diff erent expressions. Obviously, the lin- 
ear expression (F2| is sufficient for the a-edge state while 
dFl 



shows significant deviations from the numerical re- 
sults for the /3-edge state. 

The spin splitting Aeso(k) calculated from the effec- 
tive model is in better agreement with the numerical re- 
sults than the spin-orbit strengths Tg(k) (see Fig. |3j, 
even for the linear model of the /3-edge state, because 
this quantity is an average over the intrinsic and Rashba 
term of the spin-orbit Hamiltonian. 

For infinitely wide ribbons (W — * 00), the normal- 

(F3) ization factors vanish at K,K'. Thus, the spin- 

orbit splitting at K,K' is not enhanced by the Hamil- 
tonians ( F2|F4 |, compared to pure graphene. In finite 
size ribbons, however, the normalization of the edge state 
wave functions is significantly different from A/T if 
£fe > W. The edge state of a finite size a/3-ribbon 
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FIG. 9: Spin-orbit strength T e (k), 9 = from the ef- 
fective model (dashed lines) and from our extended tight- 
binding model (solid lines). The green (blue) line represents 
the Rashba (intrinsic) term of the spin-orbit Hamiltonian and 
the nearby dashed lines represent the corresponding terms 
from the effective models. The effective Hamiltonian of the 
a-edge state is Eq. |F2| . The effective Hamiltonian of the 
/3-edge state on the left-hand side (the regi on l abeled by lin- 
ear fit) is given by the linear version Eq. ( F3 1 while on the 



right-hand side (the region labeled by 5th order fit), Eq. (F4| 
is shown. Part (a) shows a calculation for a flat interface 
(zj — 0) and part (b) shows a calculation for a maximally 
tilted interface (zj = zq). 



with N unit cells in the y-direction can be written as 

|V>o(fc)) = 



i-KI 5 



i-M 2(Ar+1) 

N , v 

E ex P \~^+ m( t ) )<*,b|0>- ( F5 ) 

n=0 V ' 



At K,K' the absolute value of this wave function is con- 
stant in y-direction. Thus, the amplitude at the interface 

is not proportional to i^My!^ = in this finite-size 
case, but rather proportional to 



i-KI 2 



k^K.K' 1 

* N + l' 



As a result, the GG interface-induced SOI near K,K' is 
proportional to W~ x . Thus, the SOI of bulk graphene 
becomes dominant at K,K' in the limit W — * oo. 

In Fig. [TO] we compare the spin-orbit splitting which 
is due to the GG interface, calculated from the effective 
models ( F2 1 and ( F3 1 , with the conventional spin-orbit 
splitting which is due to the bulk graphene SOI. Obvi- 
ously, the magnitude of the interface-induced spin split- 
ting is much larger than the conventional spin splitting 
for k y^K.K'. However, at K,K' the interface-induced spin 
splitting vanishes because the edge states become delocal- 
ized over the whole graphene region so that the interface 
contribution to the spin splitting vanishes. This means 
that the spin splitting due to the bulk graphene SOI is 
dominant at K,K'. 
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FIG. 10: Comparison between the interface-induced spin- 
orbit splitting calculated from the effective model (dashed 
red line) and the bulk graphene induced spin-orbit splitting 
(green line) the edge states. The bulk graphene induced spin- 
orbit splitting is scaled up by a factor of 100. 



This is important for the QSVHE in that the Fermi 
energy must be tuned exactly to the spin-orbit gap which 
is due to the bulk graphene SOI. Qualitatively, this bulk 
SOI can be included in the effective model by adding 
the bulk graphene- induced spin-orbit splitting of the edge 
states at KK^ 



H 



SO,a/(3 



H 



eff 

SO,a/(3 ' 



A a z 



(F7) 



where 2A — 1/^eV is the spin-orbit gap at K,K' in bulk 
graphene and r 3 = ±1 for K and K', respectively. The 
red, dashed curve in Fig. [10] calculated with the ex- 
tended Hamiltonian (F7) would differ in that it is not 



(F6) exactly zero at K,K' but 2A C 
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